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The Pressure Distribution on the Head of a Shell Moving at 

High Velocities. 

By L. Bairstow, F.B.S., R H. Fowler, and D. E. Hartree. 

(Received December 12, 1919.) 

§ 1. Introduction. — The pressure at any point of a body moving through a 
fluid may be considered to consist of two components — the static pressure, 
which would be the pressure at that point, if there were no relative motion 
between the body and the fluid — and the dynamic pressure due to the relative 
motion. The sum may be called the total pressure at the point, and is, of 
course, always positive ; the dynamic pressure alone may be negative. 

Experiments, of an introductory kind, were undertaken by the authors, at 
the request of the Ordnance Committee, to attempt to determine, as a 
function of the velocity of the shell, the dynamic pressure at various points 
on the head of a shell, whose axis of symmetry and direction of motion 
coincide or nearly coincide. The results appear to be not without general 
scientific interest, and are therefore presented in this paper by permission of 
the Ordnance Committee. The authors' best thanks are due to the officers of 
the Experimental Department, H.M.S. " Excellent," on whose range the 
firing trial was carried out ; to the other officers of the Munitions Inventions 
Department, who shared the work of observation and computation ; to the 
Director of the National Physical Laboratory, where the comparative 
observations at low velocities were made in a wind channel; and, above 
all, to Mr. W. J. Goudie, of University College, London, on whose determina- 
tions of the rates of burning of powder- train time-fuzes at high rotational 
speeds the whole experiment was based. 

The fundamental idea of the experiments is the use of a service time-fuze as 
a manometer to determine the pressure under which the powder is burning.* 

The rate of burning of a service time-fuze is known from laboratory 
experiments to be affected by the total external pressure on the vents. Let 
us assume for the moment that this is the only factor that influences the 
rate of burning. A series of observations, for shells fired along the same 

* A service time-fuze contains a train of gunpowder, which is ignited by a detonator 
pellet on the shock of discharge from the gun. The " setting'' of the fuze can be varied 
so that a length of powder train depending on the setting is burnt before the magazine 
of the fuze is ignited and the shell exploded. The setting is specified by a number 
which defines the length of composition burnt on an arbitrary scale. The time of 
burning is taken as equal to the time interval between the tiring of the gun and the 
bursting of the shell. 
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trajectory at short intervals of time, determines a relation between time of 
burning and length of fuze burnt. By numerical differentiation of the 
observations, we can obtain a relation between the rate of burning and the 
time, and therefore, by comparison with the laboratory experiments, between the' 
external pressure and the time. If the circumstances of the shell's motion are 
known, this leads at once to a relation between the dynamic pressure on the 
fuze vents and the time, and so to the required relation between the dynamic 
pressure, the velocity of the shell, and the prevailing atmospheric conditions. 

The external total pressure is not the only factor that affects the rate of 
burning. Under identical physical conditions, the rate of burning varies 
from batch to batch of powder ; it is largely affected by the spin of the 
shell, and is sometimes appreciably dependent on the temperature of the fuze 
mixture. In spite of these effects, however, it became increasingly clear, a& 
experience of time-fuzes increased, that, with suitably designed control 
experiments to determine the rate of burning of the fuze in the proper 
physical conditions, a time-fuze could be used as a manometer to determine 
rough quantitative values of the external pressure during the flight of a 
shell. 

§ 2. General Scheme of an Experiment Based on the Fuze Manometer. — In the 
application described in this paper, shells were fitted with caps entirely 
enclosing the fuze ; each cap had a series of holes drilled in it at the same 
distance from the nose (fig. 1). The external pressure on the vent holes of 
the fuze is the pressure inside the cap; the latter pressure is practically 
equal to that immediately outside the hole (see § 8). Proceeding as described 
in § 1, we can obtain a relation between the pressure inside the cap and the 
velocity of the shell along a given trajectory, for a given position of the ring 
of holes. 

On the assumption that the presence of a hole has an unimportant effect,* 
we can determine the pressure-velocity relation for a series of points on the 
head of a shell by using caps with rings of holes at a series of positions, and 
draw curves (fig. 2) showing the variation with velocity of the dynamic 
pressure at a given point. For any given velocity we can draw curves 
(fig. 3) showing the distribution of dynamic pressure over the head of the 
shell. By integrating these dynamic pressures, we can deduce the contribu- 
tion to the drag due to the dynamic pressure on the head of the shell. The 

* In order to make use of these relations we must assume that the total pressure on the 
outside of the {small) hole is the same as the total pressure at the same point of the shell in the 
absence of any hole. This is the assumption made in all aerodynamical experiments on 
pressure distribution at low velocities. Its application to high velocities is at present 
incapable of justification, either experimental or theoretical. We shall return to this 
point in § 8. 

Q 2 
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total drag must be made up of the contributions of (1) the dynamic pressure 
on the head, (2) the (negative) dynamic pressure on the base, and (3) skin, 
friction. 

If we could simultaneously determine the total drag and the two contribu- 
tions of the dynamic pressure over a range of velocities, we could probably 
make some advance to an understanding of the mechanism of the resistance 
of a compressible fluid at high velocities. 

It seems probable, from general considerations, that the pressure on the 
exact nose of a shell will obey closely the formula given by Lord Kayleigh.* 
Caps with a hole at the nose precisely were not used in this trial; they 
could, however, be so used, and a direct test of Bayleigh's formula could be 
made, which would be of great interest. The curves shown in fig. 2 are of." 
the same general shape as Kayleigh's curve for the point of the nose, but the^ 
pressure falls off so rapidly near the nose that a precise comparison on: 
present results is not possible. 

A number of difficulties, many of which will occur to the reader, have 
been tacitly ignored in the foregoing sections. These will be discussed later 
on in § 8. 

§ 3. Details of Firing Trials. — An 18-pdr. field gun (on a high-angle- 
mounting) was selected for the trial ; the muzzle velocity, determined^ 
during the experiment, was 1603 f.s. The shells used are shown in fig. 1.. 
The caps were rounded off at the apex, and the term " virtual point " will- 
be used to denote the position of the point of a cap with a perfectly sharp- 
point. 

There were five series of caps, each with a ring of six holes. Four series- 
had holes of diameter 0*05 inch and one series 0*1 inch; these series will be 
referred to by letters as follows : — 

Table I. 



Series letter. 


Number of 


Diameter of 


Distance of holes from virtual point of cap, 


rounds. 


holes. 


measured along the axis of the shell. 






inches. 


inches. 


A 


30 


0*05 


1-30 


B 


50 


0'05 


2-55 


C 


30 


0-05 


3-80 


c* 


30 


o-i 


3-80 


D 


30 


0'05 


4'80 



Series C* was inserted as a control, to determine whether the effect of the 
size of the hole on the pressure was likely to be considerable. 

* " Aerial Plane Waves of Finite Amplitude," l Scientific Papers,' vol. 5, p. 610. 
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Groups of five rounds were fired at fuze settings 3, 6, 9, 12, 16, 20 at an 
elevation of 20° in each series. The extra rounds in series B were fired at 
50°. The majority of the rounds in series B were fired by day, and the 
position of the shell bursts was observed. The remaining rounds of all series 
were fired by night. Times of burning of the fuzes were recorded for all 
rounds fired ; the means for each group are given in Table II. 

Table II.-— Mean Observed Time of Burning for each Group, and Average 
Kesiduals* of the Eounds in each Group, with the Becalculated Times 
of Burning. Times in seconds. Zero point taken as +0*24 fuze 
settings. (See § 4.) 



Fuze 
setting. 


Series A, 20°. 


Series B, 20°. 


Series B, 50°. 

■ 


Recalcu- 
lated. 


Ob- 
served. 


Average 
residual. 


Becalcu- 
lated. 


Ob- 
served. 


Average 
residual. 


Recalcu- 
lated. 


Ob- 
served. 


Average 
residual. 


3 

6 

9 

12 

16 

20 


2*41 

5'15 

8*04 

11-07 

15-26 

19-47 


2-47 

5-28 

8*02 

10*90 

15*26 

19-50 


0-15 
0-08 
0*13 
0-21 
0*12 
0-12 


2-48 

5-02 

8-32 

11-48 

15-78 

20-07 


2-42 

5*34 

8-36 

11-47 

15-82 

20-16 


0*04 
0-12 
0'06 
0-35 
0-15 
0*12 


9*03 
12-85 
18-23 

23*82 


9*01 
12-73 
18-19 
23-68 


0*16 
0-08 
0-16 
0-27 




Series C, 20°. 


Series 0* 20°. 


Series D, 50°. 


3 

6 

9 

12 

16 

20 


2-67 

5-68 

8*81 

12-04 

16-45 

20-84 


2-70 

5-65 

8-92 

12-02 

16-46 

20-87 


0-07 
0*12 
0-29 
0-18 
0-34 
0*29 


2-68 

5-79 

9-14 

12-57 

17 -09 

21-53 


2-70 

5-72 

9-12 

12-57 

17-15 

21-38 


0-04 
0-15 
0-14 
0-08 
0-10 
0-31 


2-84 

6-02 

9*46 

13*02 

17*73 

22-28 


2-86 

5*97 

9-51 

13*08 

17-64 

22-31 


0-09 
0-17 
0-32 
0-22 
0-36 
0-17 



The following abnormal rounds are omitted from means and average 
residuals (see § 8) : — 

Series A: Set 20, 10-10 ; Set 16, 7*03; Set 12, 8-32. 
Series B : 20°, Set 12, 5*7 ; Set 3, 1*92. 

50°, Set 16, 16-89 ; Set 12, 8*03. 
Series C : Set 12, 10*23 ; Set 9, 4*62, 7*81 ; Set 6, 3*98. 
Series D : Set 6, 2-83 ; Set 3, 1-84. 

The times were recorded to the nearest 0*01 second, but the time for any 
one round is not accurate to this degree. Experience has shown that under 
favourable conditions, such as obtained in this experiment, and with three or 

* The arithmetic average difference of the members of the group from the arithmetic 
mean of the group. 
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more observers, the time for any round can be determined with a probable 
error of 0'04 second. 

By means of the observations of position, the ballistics of the shells were 
determined. The conditions in the ripper air were known from aeroplane 
observations. The wind structure was also observed. From these observa- 
tions and routine calculations which need not be specified, the following 
quantities were determined and tabulated for the position of the shell at 
every integer second after gunfire along both the 20° and 50° trajectories :— 

(1) The relative velocity of the shell and the air, v. 

(2) The static pressure of the air, p. 

(3) The density of the air, p. 

(4) The ratio v/a, where a is the velocity of sound in the air. 

§ 4. Laboratory Determination of the Bate of Burning of the Fuzes. — In order 
to determine the rate of burning of the fuzes used under the proper physical 
conditions, it is essential that laboratory experiments should be carried out 
with fuzes filled with the same batch of powder and under the same filling 
conditions as the fuzes used in the gun trials. Euzes fulfilling these con- 
ditions may be described as fuzes of the same lot. The rate of burning of a 
fuze of given lot at given pressure and temperature varies largely with the 
spin, especially at low pressures. It is therefore essential to carry out the 
laboratory experiments with the fuze spinning at the correct rate, and (some- 
times) at the correct temperature. As the velocity of the shell cannot affect 
the fuze, and the deceleration is small compared to the radial acceleration in 
the powder channels, the physical conditions of an actual firing trial are then 
properly represented in the laboratory experiment. 

The initial spin of the shells in the trial was 11,660 revolutions per 
minute, and it may be assumed that it remains roughly constant for the 
fairly short times of flight with which we are concerned. In Table III the 
rates of burning of the fuzes used are given, determined at a spin of 
11,750 revolutions per minute. This determination was effected by Mr. 
Goudie by coupling a dummy fuze body, holding the rings containing the 
powder train, to the shaft of a de Laval turbine ; the dummy was run in an 
air-tight box with an automatic control of the pressure. 

Experiments were also carried out to determine the effect of temperature 
and change of spin. Both these effects were negligible for the present trial. 

The arbitrary graduations defining the fuze setting on these fuzes run 
from to 22, and the length of powder burnt varies linearly with the fuze 
setting on this arbitrary scale. Consequently, under constant external 
conditions, the time of burning is a linear function of the graduations or fuze 
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Table III. — Showing the Eelation between the Eate of Burning and the 
Pressure on the Vent Holes for the Fuzes used in this Experiment, 
determined by Mr. W. J. Goudie at University College, London. 
Spin. — 11,750 revolutions per minute. 

The rate of burning is given in fuze settings per second. The pressure is 
given in atmospheres. 



Eate of burning. 


Pressure. 


Difference. 


Rate of burning. 


Pressure. 


Difference. 


0-5 

0-55 

0-6 

0-65 

0-7 

0*75 

0*8 

0-85 

0-9 

0-95 

1-0 

1-05 

1-1 

1-15 

1-2 


0*407 
0-464 
6 -522 
0-581 
0-641 

0-703 
0-767 
0-833 
9-901 
0-971 

1*044 
1-120 
1-200 
1*285 
1-376 


0-057 
0*058 
0-059 
0-060 
0-062 

0*064 
0-066 
0*068 
0-070 
0-073 

0-076 
0-080 
0-085 
0-091 
0-096 


i 

1-25 

1*3 

1-38 

1-4 

1-45 

1 1-5 
1 1 55 
1-6 
1-65 
1-7 

t 

j 1-75 

i 

i 

l 
i 


1-581 
1-699 
1*829 
1-974 

2*137 
2 *o2JL 
2 *528 
2*763 
3*030 

3*335 


0*107 
0-118 
0-130 
0-145 
0*163 

0-184 
0*207 
0-235 
0*267 
0*305 



setting except for very short times. The rates of burning of Table III are 
given in fuze settings per second, and apply to all stages of the combustion 
under constant external conditions. In order to construct from the obser- 
vations, as accurately as possible, a curve of fuze setting against time of 
burning (see fig. 4, § 8), it is important to determine the constant in the above 
linear relation. This constant may be regarded as the fuze setting which 
burns for zero time and is called the zero point of the fuze. It is not 
necessarily at the zero graduation. Its position depends on the internal 
arrangement of the powder-train connections, and it must be determined 
by experiments on the actual fuze bodies used. It appears that the zero 
point of the fuzes used was +0*24 settings. It is assumed that this zero 
point remains unaltered when a fuze is fired from the gun. If this assump- 
tion is incorrect, the rates of burning observed during firing trials at short 
times will be in error, but those at later times will be unaffected. 

§ 5. Analysis of the Observations.— -The data provided by the observations 
described in §§3, 4, are sufficient, on analysis, to provide information as to 
the value of the dynamic pressure for values of v/a between 0*7 and 1*4* for 

•* Note added, March 10. — Owing to the effect of an uncertainty in the zero point 
(see § 8), the values of the dynamic pressures for values of v/a greater than 1*2 must 
be given less weight than those for values of 1*2 and less. 
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each position of the holes in the cap. It is not necessary to describe in 
detail the numerical methods used. After the final results for the dynamic 
pressures were obtained, the inverse process was gone through of integrating 
the corresponding rates of burning to obtain the time for each fuze setting 
us a check on the accuracy of the work. These recalculated times are given 
in Table II, and it will be seen that the agreement with the observed time 
is satisfactory for all series, the difference between the two for any one group 
being less than half the average residual for the observed group in nearly 
all cases. 

In series B, which was fired at 20° and 50°, the recalculated times of 
burning are systematically slightly too short at 50° and too long at 20° at 
long settings. The differences, however, are not bigger than those for other 
series. 

In the usual aerodynamical notation, the dynamic pressure may be put in 

the form 

pv% (v/a), 

where v is the velocity, p is the air density, a is the velocity of sound, 
and f p is a numerical coefficient known as the pressure coefficient. This 
coefficient f p may (by the theory of dimensions) depend on other arguments 
besides v/a, but it is unlikely that any others are of serious importance in 
this connection. It is known that f p is roughly constant at low velocities ; 
this eonstancy probably persists nearly up to v/a = 0*7. If pressure, 
density, and velocity are measured in any self-consistent set of units, f p is 
independent of the set of units chosen, and is more suitable than the 
dynamic pressure itself for tabulation as it does not depend on p. The 
values of f pt determined by these trials were calculated from the dynamic 
pressures, and are tabulated in Table IV as functions of v/a. They are 
-also shown graphically in fig. 2. 

Values of the pressure coefficients, determined in a wind channel at 
the National Physical Laboratory, at a velocity of 40 f.s., are also given 
in Table IV for the same positions of the holes along the axis of the shell. 
The dynamic pressures at this velocity were determined more exhaustively 
than this along the whole shell. It will be seen in Table IV or fig. 2, 
that the observed values of f p (0*7) fit on satisfactorily to the wind channel 
values. 

By plotting the values of f p in Table IV for given values of v/a, opposite 
the position of the holes at which f p was determined, a diagram is obtained, 
fig. 3, showing the distribution of dynamic pressure (coefficients) along the 
shell at high velocities. The curves are drawn for values of v/a equal to 
1-47, 1*2, 1-0, 0-9, 0-75, and 0*04 (wind channel velocity). The values of the 
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•8 



' 7 A 

y/a^l'2 

•6 

V/a*075 



Fig. 3. 
Showing Distribution of Pressure along the 18-pdr. Shell 

of fig. 1 at Various Velocities. 

Continuous curves ( ) are drawn through points, 

marked according to the values of vja to which they refer, 
obtained from firing trials, and, for the nose of the shell, 
from calculation. 

Dotted curves ( ) are extrapolations of the 

continuous curves, which are needed for integrating up the 
pressure on the head. 

Dash and dot curve (— . — . — . — ) is drawn through 
points marked x , determined by wind channel experiments, 
at a velocity of 40 f.s., or value of (vja) of 0*04. 

The section of the shell below js drawn to half size, and 
the longitudinal scales of the two diagrams are the same. 
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pressure at the nose, at which these curves start, and the curve for the 
pressure at the nose, given in fig. 2, are calculated by Kayleigh's formulae 
<(loc. cit). 
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In fig. 3 the exact shape of the curve between the nose and the position 
of the holes in series A is naturally not very definite. This is not important 
for the purpose of § 6, but it could be wished that the curves were determined 
by observation behind the position of the holes in series D. 

§ 6. Integration of the Pressure over the Head of the Shell, — It can easily be 
shown that if p (y) is the dynamic pressure at the point on the shell at which 
the distance from the axis is y, then the total force due to the dynamic 
pressure on the shell in its symmetrical position is along the axis and is 
given by the expression 
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y p (y) <ty> 







where r is the radius of the shell. The mean component of the dynamic 
pressure on the head, parallel to the axis of the shell, is therefore 
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and the coefficient of this contribution to the drag, or mean dynamic 
pressure coefficient f p (v/a), satisfies the equation 

This coefficient has been calculated for several values of v/a, and the results 
are given in Table IV and plotted in fig. 2. 

Again, if we express the total drag E on the shell in the form 

E = Apv 2 f n (v/a), (2) 

where A is the area of the cross-section of the shell, then the drag per unit 

area is given by 

E/ A = pv% (v/a), 

The coefficient f R (v/a) is directly comparable to f p (v/a), the coefficient of 
the mean dynamic pressure over the head, which latter is one of its main 
contributors. 

The values of / R , the coefficient of the drag per unit area, are not known 
really accurately for shells of this shape. They certainly lie, however, 
between the two sets of numbers given in Table IV and the two curves of 
fig. 2, approximating to the lower values at high velocities and the higher 
values at low velocities. 

The pressure coefficient f B (v/a) } corresponding to a total vacuum over the 
base, is clearly given by the equation 

/b (v/a) = p/(f*P) = 7" 1 (v/a)", (3) 

where 7 is the ratio of the specific heats of the gas. Values of f^(v/a) are 
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also tabulated in Table IV and plotted in fig. 2. It will be seen that the 
coefficient / K is less at all velocities here reached than the coefficient f B alone. 
The vacuum on the base is therefore by no means complete at a velocity 
represented by v/a = 1*5. 

It would be of very great interest if the precise degree of vacuum could be 
properly determined by experiment. This would not be possible by the use 
of a powder train fuze, for no such fuze will burn at all at the necessary low 
pressures. 

§ 7. An immediate value attaches to the results now recorded, as they 
provide an adequate explanation of the complex behaviour of fuzes. The 
changes in rate of burning which accompany change of muzzle velocity 
or change of diameter of shell are all explained qualitatively by curves of 
the kind shown in fig. 2. It is now possible to predict the behaviour of fuzes 
before trial in a manner, and with a degree of certainty, unattainable prior to 
the war. 

It is interesting to discuss briefly the records of pressure-distribution on 
the head of a shell in connection with other knowledge bearing on the point. 
The curves of fig. 2 show that at all speeds in the range covered the pressure 
has a maximum positive value at the nose of the shell. The pressure falls 
rapidly as the point of observation moves towards the base, and is negative 
some distance before the cylindrical part of the shell is reached. The deter- 
mination of the pressure at the nose has only been made at low speeds 
(v/a = 0*04), but is then found to be independent of the viscosity of the fluid 
medium. At these low velocities, calculations and experiments have shown 
that this independence of viscosity extends so far as sensibly to cover all 
positive dynamic pressures. The maximum observed pressure at the nose 
agrees with that calculated on the basis of an in viscid fluid, and the degrada- 
tion of energy caused by viscosity shows an increasing effect as the point of 
observation recedes from the nose. The complete recovery of pressure on the 
base of a shell which is indicated by the theory of an inviscid fluid is not even 
approximately represented by the experimental value. 

The reduction of the observations of the trial shoots has been made on the 
hypothesis that the effects of viscosity on the pressure on the head of a shell 
are negligible, and the results are quite consistent with that hypothesis. It 
may be expected, therefore, that the pressure could be calculated approxi- 
mately from the equations of motion of a body in an inviscid compressible 
fluid. In the event of a solution of those equations, the observations now 
recorded would provide a suitable check on the soundness of the initial 
assumptions. 

Viewed as an addition to the experimental investigation of the motion of 
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shells, the trials have an importance which is too obvious to require detailed 
notice. 

§ 8. Possible Disturbing Factors, and Accuracy of the Eesults. — In conclusion 
it may be well to return to a number of difficulties which have been more or 
less ignored in the earlier sections. We have stated that the pressure at the 
outside of the holes in the cap is practically the same as the pressure on the 
vent-holes of the fuze inside. The factors that might be thought to destroy 
this equality are (1) the centrifugal force in the rotating column of air 
between the vent-holes and the outside of the cap ; (2) the pressure difference 
necessary to exhaust the products of combustion through the holes at the 
proper rate ; (3) an initial delay while the interior of the cap is filling up to 
the true pressure outside. The effect of (1) is easy to calculate. In the case 
of (2) the maximum rate of efflux required is known from the composition of 
the powder and rate at which it is burnt. Observations were made of the 
pressure differences required to attain a measured rate of efflux with the 
actual caps. These observations also apply to (3). As a result it appears that 
the total pressure difference will only amount to more than 0*01 atmospheres 
for a very short time (of the order of 1/20 sec.) at the start. The effect is 
therefore completely negligible. 

A far more serious question is, whether the pressure on the holes is actually 
the same as the pressure on the surface of the cap in the absence of the 
holes, the holes themselves being small. This is an accepted principle in 
aerodynamical work at low velocities,* and there seems no valid reason to 
expect it to fail at higher velocities. As a control of this assumption, the 
comparative series C and C* were fired ; in C* the diameters of the holes 
were twice those of C. There is definite evidence that the dynamic pressures 
determined as above are different for these two series. The difference is not 
appreciable at the higher velocities, as the times of burning at setting (3) are 
the same, but for values of vja between 1*1 and 0*8 the pressure determined by 
C* is distinctly less than the pressure determined by C (see fig. 2). The holes 
of the C* series are therefore not " small holes " at these velocities. It thus 
remains somewhat doubtful whether the holes of series A, B, C, D are really 
" small." If they are not, then the true dynamic pressures are probably 
somewhat higher. The difference is likely to be small, but it is not possible 
to assign to it a numerical upper limit. 

* The best direct experimental evidence will be found in the following papers : 

(1) G. Fuhrmann, " Theoretische und experimented Untersuchungen an Ballon - 
modellen," 'Jahrbuch der Motorluftsehiff-StudiengesellschafV vol. 5, p. 12 (1911); 

(2) B. M. Jones and C. J. Paterson, " Investigation of the Distribution of Pressure 
over the Entire Surface of an Aerofoil," Advisory Committee for Aeronautics. 
Reports and Memoranda, No. 73, March, 1913. 
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The remaining points are less important. We have tacitly assumed that, 
throughout the motion, the axis of symmetry of the shell and the direction of 
motion of its centre of gravity coincide. The departures from this position of 
symmetry will be very small. They may exceed 1° or 2° initially for a short 
time, but for most of the motion they will be much less. The effect of such 
departures on the average pressure inside the cap has been determined at low 
velocities and is known to be very small indeed. The whole effect will 
prooably be negligible here. 

Finally, we come to the uncertainties in the behaviour of the fuze, and 
especially in the position of the zero point. All the fuzes used were subjected, 
as far as possible, to precisely similar treatment. There is no room for any 
suspicion of systematic differences between different parts of the material 
used. Eandom differences alone should remain, and the extent of these is 
shown by the average residuals of Table II. It will be observed also that a 
number of rounds are , omitted from the means and mean differences of this 
Table. Since all these omitted rounds are short, and since they occur in all 
series except C* with its large holes, some occasional mechanical trouble, such 
as the blocking of the holes by slag or paint, is suggested. 

On the other hand, an error in the zero point will cause systematic errors 
in all the dynamic pressures at the higher velocities. This is shown by fig. 4, 
where it is obvious that the precise position of Z will seriously affect the slope 
of the time-fuze setting curve for short times. It is by determining the 
slopes of such a curve that the rates of burning are obtained. 




= Observations^ of 
time of burning 
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The final results have all been given to a nominal accuracy of 0'005 inf p . 
This corresponds to an actual pressure of about 1/200 of an atmosphere at; 
vol. xcvii. — A. R 



218 Captain J. W. Bispham. Distribution of the 

the velocity of sound. It must not be supposed that this represents the real 
accuracy of the results, especially at velocity ratios above 1*1 or 1*2, where 
the effect of the zero point may be serious. It does seem, however, possible 
to claim that the tables and curves of this paper give values of these coeffi- 
cients of the true order of magnitude, and show correctly the type of change 
that takes place as the velocity of the shell passes through the velocity of 
sound. 



An Experimental Determination of the Distribution of the 
Partial Correlation Coefficient in Samples of Thirty. 

By Captain J. W. Bispham, E.E. 
(Communicated by Prof, C. J Martin, F.RS. Eeceived November 26, 1919.) 

Part I. — Samples from an Uncorrelatecl Universe. 

1. Introductory: 

In 1908* " Student " dealt experimentally with the distribution of the 
total correlation coefficient of small samples. In particular, he dealt with 
values of n as low as 4 for the case of zero correlation in the sampled 
population. In 1913 H. E. Soperf theoretically determined the mean correla- 
tion and the standard deviation of the distribution of correlations to second 
approximations. In 1915 E. A. Fisher J gave an equation for the frequency 
distribution of r, and in 1917, as a result of a co-operative study by 
H. E. Soper, A. W. Young, B. M. Cave, A. Lee and K. Pearson^ this was 
reduced to suitable form for numerical manipulation, and the frequency 
distributions and frequency constants for samples of size ranging from n = 3 
to n = 400 were given for values of the correlation in the sampled population 
ranging from p = to p = 0*9. The present experimental investigation was 
commenced in 1914, but had to be put aside during the war. It was intended 
to determine whether the distribution of partial correlation coefficients for 
samples as small as 30 showed greater dispersion than is observed for total 
correlation coefficients. Yule|| has shown that for normal distributions and 
large samples the standard deviations of the distributions should be of the 

* ' Biometrika,' vol. 6, p. 302, et seq. 

t ' Biometrika/ vol. 9, p. 91, et seq. 

\ l Biometrika,' vol. 10, p. 507, et seq. 

§ ' Biometrika,' vol. 11, p. 328, et seq. 

|| G. U. Yule, 'Roy. Soc! Proc./ A, vol. 79 (1907). 



